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What has become the standard approach to relativisitic Coulomb excitation (RCE) was 
| proposed by A. Winther and K. Alder (WA) in 1979 The projectile nucleus is assumed to 
travel along a straight-line orbit parallel to the z axis, with impact parameter b, at constant 



speed v. The magnitude of the impact parameter is large enough so that nuclear interactions 
between the target and projectile are negligible. Because of the assumed large projectile 
momentum, the electromagnetic impulse the projectile receives due to its interaction with 
the target has little effect on its trajectory, so the projectile maintains its constant speed and 
impact parameter throughout the collision. As the projectile passes, the target nucleus feels 
the time-dependent projectile electromagnetic fields, which induce transitions between the 
quantum states of the target. Starting from the assumption that the target is in its ground 
state at t — — oo, WA used first-order perturbation theory to calculate the occupation 
probabilites of excited target states at t = +oo, and used these probabilities to obtain 



Coulomb excitation cross-sections 



21] 
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An important ingredient in the calculation of target transition probabilites is the inter- 
action potential felt by the target as the projectile moves past it. WA took this to be the 
classical electromagnetic field of the moving projectile. A justification of this assumption 
can be found in the work of Alder, Bohr, Huus, Mottelson and Winther (ABHMW) 
These authors used the lowest-order of perturbation theory to calculate the target transi- 
tion amplitude as a result of photon exchange with the projectile, in a situation in which 
the projectile motion was described by quantum mechanics. They found that this photon- 
induced transition amplitude was the same as the target transition amplitude induced by 
the classical electromagnetic field of the projectile, again calculated in the lowest order of 
perturbation theory. 

It is well known |5j that the quantum-mechanical treatment of the interaction of two 
charged particles yields the same result whether the interaction between them is 

• The classical electromagnetic field, or 

• The quantized electromagnetic field, calculated up to terms of order e 2 . 

This is essentially the same result as reported by ABHMW0]. 

In the 25 years since the WA paper, attempts have been made to improve their calculation 
of transition amplitudes by going beyond first-order perturbation theory ^ , 12 , , 14, 
Q Q, 0, Q- In particular, the time- dependent Schrodinger equation that governs the 
occupation amplitudes of the target states has been integrated numerically from t = — oo to 
t = oo as a set of coupled-channel equations, within a finite set of target basis states. In these 
calculations, the classical electromagnetic field has been used as the interaction potential. 
But this clearly goes beyond the justification provided by the work of ABHMW, which 
established the connection between the classical and quantized field only up to first-order 
perturbation theory. 

In this paper we consider a somewhat restricted problem. We will assume that transition 
charge and current densities that characterize the projectile are specified function of position 
and time, pp(r, t), Jp(r, t). In fact this is what is usually done in RCE calculations, with 
pp(r, t) and Jp(r, t) taken to be the charge and current densities appropriate to a spherically- 
symmetric charge moving with constant speed v along a straight-line trajectory with impact 
parameter b. We will not restrict the trajectory or the shape of the projectile; we only 
require that pp(r, t), Jp(r, t) be specified functions j^j]. This means that we are neglecting 
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the effects of projectile-state changes on the coupling of the projectile to the electromagnetic 
field. We will see that this commonly made assumption allows us to make a much stronger 
statement about the connection between the classical and quantized-field treatments of the 
excitation of the target. In fact, we will prove that they yield the same results to all orders 
of perturbation theory, not just to first order. 

In Section II we will develop the basic expression for the transition amplitude in terms 
of the time-dependent interaction potential. In Section III we will evaluate this expression 
in the particular case in which the interaction is provided by the quantized electromagnetic 
field, calculated in the Coulomb gauge, and in Section IV we will show that this result is 
precisely the same as if we had used a classical electromagnetic field for the interaction. In 
Section V we will show that this agreement also holds if we had used the Lorentz gauge 
for both the classical and quantized field calculations. A more accurate treatment of the 
coupling between the projectile and the electromagnetic field is discussed in Section VI. 



II. THE TRANSITION AMPLITUDE 

Consider the time-dependent Schrodinger equation 

^ = [ho + v(t) } $ (II.l) 

where ho does not depend explicitly on t. This equation is expressed in the interaction 
representation by expanding tp in terms of the normalized eigenstates of h : 

^(t)=^e- i ^ t a 7 (t)$ 7 , (II.2) 

7 

with 

h % = £ 7 $ 7 (II.3a) 

< $ 7 | > = 5^ (II.3b) 

If $ a is the initial state, which existed when t — > — oo, then 

a 7 (-oo) = 5 7>a . (II.4) 

The transition probability to state $ 7 is given by |a 7 (+oo)| 2 . 
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Because of equation the a a (t) obey the set of coupled differential equations: 

iM^t) = e *° 7/3 * W)\if> a /?(*)> ( IL5a ) 



with 



f2 7/3 is defined by 



= <%wm P >. (ii.5b) 



tor/, = (H.6) 



The integral equation equivalent to equation (|II.5aj) . incorporating the initial condition 
(flL4|) . is 

/* dt' 
-OO 1 < L Q 

This can be iterated to develop a perturbation series in powers of v(t): 

/* rli' r f rli' r 1 ' rli" 

-co ^ ^ J —OO ^ J —OO ^ 

(II.7) 

The amplitude for a transition from <3> a at t = — oo to $ 7 at t = oo is a 7 (oo). Thus the 
perturbation series expansion of the transition amplitude can be written 

POO lA. I>t n JJ. 

a 7 (oo) = ^ a + J2J2 / ^^v„(t n ) ^e^" 1 ^^)... (II.8) 

__1 a </-00 ^ J -OO % ' 1 



n=l /3...T ' 



oo ^ J-oo in J-oo in 

III. SPECIALIZATION TO THE QUANTUM THEORY OF RELATIVISTIC 

COULOMB EXCITATION 

In the usual approach to relativistic Coulomb excitation (RCE), the target nucleus is at 
rest in the reference frame. The rapidly moving projectile nucleus follows a prescribed orbit, 
which we will leave unspecified. The projectile interacts with the electromagnetic field via 
its time-dependent charge and current densties pp(r,t) and Jp(r, t). 

The Hamiltonian ho of Equation (jll.l|) will refer to the target internal degrees of freedom, 
plus the free electromagnetic field. An eigenstate of ho will therefore be specified by a target 
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state 0/3, plus a specification of the numbers of photons in each of the quantized modes of 
the field. If the energy of the target state 4>p is ep, then a transition from target state 4> a to 
(f)p with creation of a photon of momentum ftq will result in an energy increase of 

{ep + hcq) - e a , 



so that the quantity Qp a referred to in Equation ()H.8j) is 

(e p + hcq) -e a e p -e a 
llp a = = h cq = Up a + cq. 

Similarly, if the <fr a —>■ <pp transition were accompanied by the absorption of a photon of 
momentum ftq, we would have 

(ep - hcq) -e a ep - e a 
ilp a = = — cq = ujp a - cq. 

Since the calculation involves the electromagnetic potentials, we must choose a gauge. 
In this section and the next, we use the Coulomb, or radiation, gauge, j^, 0, Q] because it 
allows the simplest treatment of the quantized electromagnetic field. In Section V we will 
discuss the modifications required for the Lorentz gauge. In the Coulomb gauge, the vector 
potential is a solenoidal field 

V-A(r,t) = 0. (III.l) 

Here r labels points of space relative to an origin at the target center, which we assume to 
be fixed. We describe the photon field in terms of normal modes, each labelled by a wave 
vector q. Associated with each q are two polarization vectors e q i and e q! 2- The three vectors 
(q, e qi i, e q> 2) form an orthonormal coordinate system. The field operator A(r) is expressed 
in terms of photon creation and annihilation operators (a qj -,a qj ) by the expansion 

A « = c \/y EE Jh I '"'"'V".,, + e-* qr e q X, ] • (ni.2) 

q j = l V 

Here V is the quantization volume. The number of modes per unit q-space volume is V/87T 3 , 
for each polarization e qj . Since the plane waves that describe the modes are transverse (e qj - 
q = 0), the solenoidal condition (|III.1|) is automatically satisfied. With the normalization 
given in equation (jIII.2|) . the (a q -,a q j) obey the commutation relations 

[ fl q,j> a J,j' ] = <W<W (III.3a) 
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[ a q,i> fl q',i' ] = [ a t,V a q'J' 1 = (HL3b) 

Their matrix elements with respect to states with n q j photons in the mode q, j are 



< ri^j | a+j | n qJ > = ^..n^+iV^qJ + 1 (III.4a) 

< n'nj I a Q,i I n d,j > = ^.n^-iV^qJ (III.4b) 

In the Coulomb gauge, the full Hamiltonian for our system is 0, 0| 

# = h -- c J d 3 r{ J P (r,t) + J T (r) ) ■ A(r) + J d 3 rp T {r) J d 3 r' PP ^ (III.5a) 
h = Ht + H^ (III.5b) 

Ht is the target internal Hamiltonian, and 

2 

^ = J2 hc( i a ij a ^ 

q j=i 

is the free-field photon Hamiltonian. There is no term in equation (jIII.5a|) corresponding to 
the kinetic energy of relative motion. We are using the usual RCE picture of the projectile 
moving on a prescribed classical trajectory, and therefore the projectile-target relative coor- 
dinate is not one of the degrees of freedom of the problem. The projectile charge and current 
densities (pp(r',t), J P (r', £)) are to be regarded as specified functions of r',t, whereas the 
target charge and current densities are to be regarded as operator functions, to be repre- 
sented in the calculation by their matrix elements [p^^lr, t), [J] 7i/ g(r, t) with respect to the 
target eigenstates </> 7 , </>/?. 

The interaction to be used in Equation (jll.lj) is 



v (t) = -}-Jd 3 r[( J P (r,t) + J T (r) ) • A(r) ] + J d 3 rp T (r) j dV ^ ( _^ (III.6) 
It is conveniently regarded as a sum of two terms: 

v {t) = J d 3 rp T (r) J d 3 r'^^ (III.7a) 
which does not change the number of photons, and 

vi(t) = -- c J d\[ ( J P (r,t) + J T (r))-A(r)] (III.7b) 
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which changes the number of photons by ±1 because of the presence of a qj and •/ in the 
expansion of A(r,i) (Equation ()HL2|) ). 

We are interested in a situation in which neither the initial state <ft a nor the final state 7 
has any photons [23]. Thus we exclude bremsstrahlung processes. Moreover, we assume that 
every photon absorbed by the target at time t was emitted by the projectile at an earlier 
time t', and every photon emitted by the target at time t will be absorbed by the projectile 



at a later time t' 



24j . This implies that the creation and annihilation operators entering into 



the expansion (Equation (|IL8p ) for the transition amplitude will occur in pairs a, 



1,J ••■ Q q,i 



Specifically, if the t p integrand in Equation (|II.8J) contained 



e q ,j • [JtWM a qj e jqr = e qJ ■ [J T W(r) a qj e 



iqr 



corresponding to absorption of a (q, j) photon at the target at time t p , it must be that this 
photon was created at the projectile at an earlier time t' . This means that farther to the 
right in the multiple integral of Equation (|II.8J) there occurs a factor 

picqtp 

-^-^..J P (r',g a&e-** 

and t' must be integrated from — oo to t p . Similarly, if the t p integrand in Equation (|11.8|) 
contained 



i^lfj,i/tp piiu)^ -\-cq)tp 

-j- ■ [JtI^W a+,e— = - ■ [J T ]^(r) <, 



corresponding to creation of a (q, j) photon at the target at time t p , it must be that this 
photon will be annihilated at the projectile at a later time t' p . This means that farther to 
the left in the multiple integral of Equation (jll.8|) there occurs a factor 

— ^-e qii .J P (r',g a q ,,e + ^ r 

and t'p must be integrated from t p to oo. The result of these two t' p integrations is that the 
double integral over t p and t' p in Equation (jll.8|) is replaced by the single integral 

/ -JLj"^ w, u (t p ), (III.8) 
J-oo in 
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where w^ u (t p ) is defined by 



2 

47T x-^^ h 



q j=l 



+ I d 3 re qj ■ [J T ]„(r)e J (- c ^ + ^ ^ ^ | ^•J P (p' i g e < K^ 



Note that it is possible to group all the contributions to the t' p range in Equation (|II.8J) into 
these two intervals, —00 < t' < t p and 00 > t' > t p , only because of our assumption that 
the projectile current density is independent of what occurs at the projectile or target. 

We can repeat this process for every a qj . . . a^j pair ocurring in the multiple integral of 
Equation (|H.8|) . thereby replacing every t,t' integration (t ocurring at the target, t' at the 
projectile) by a single t integration. All that will remain in Equation (jll.8|) is integrals of 
the form (jlll.8|) obtained in this way, and integrals of the form 

, in 

Note that here VL^ V = lu^, since Vo(te) does not create or annihilate photons. The remaining 
combination of the integrals 

^e^w pu (t p ) and / %^v (t e ) 

is equal to 

a 7 (oo) = 5 7 , Q + E / ^ e ^ tn ^(*n) / ^e"^%(Vi) • • • (HI.10) 

00 & i-oo in in 

in which the time-dependent perturbation v(t) is defined by 

v(t) = v (t) + w(t). (ni.ii) 



Comparison with Equation pi. 8(1 shows that the transition amplitude for the excitation of 
target state 7 under the influence of the quantized electromagnetic field is precisely the 
same as if the target had experienced the time dependent effective field v(t), from which all 
photon degrees of freedom have been removed. 
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In the next Section, we will discuss the physical meaning of the effective target interaction 
v(t) defined in Equation (jlll.IIj) . 

The structure of Equation (|III.10j) can be elucidated by examining the n = 2 term in 
detail: 

|f e^ 2 [v (t 2 ) + w(t 2 )}^ / H£e*"*A [v (h) + w{t x % a (111.12) 

The ^0(^2)^0(^1) combination arises from the n = 2 term of Equation (jll.8|) . when the Vq part 
of v(t) (Equation (|III.7a|) ) is used both times. The vo(t 2 )w(ti) and w(t 2 )vo(ti) combinations 
arise from the n = 3 term of Equation (jll.8|) . when one of the three interactions is chosen 
to be Vq and the other two are chosen to be photon interactions V\ (Equation (|III.7b|l ). The 
w(t2)w(ti) combination arises from the n = 4 term of Equation (|11.8|) . with v\ acting four 
times. This involves the exchange of two photons between projectile and target, the target 
interactions occurring at t 2 , t\ and the corresponding projectile interactions occurring at 
t' 2 ,ti, with both t 2 and t[ integrated from —00 to 00. These components of the n = 2 term 
of Equation l|lll.10|l are illustrated in Figure 1,2 and 3. 



IV. THE EFFECTIVE INTERACTION 



If we define four fields <p(r, t) and A(r, t) by 
0(r,t) = I r/V 



3 p P (r',t) 



|r — r 

2 



(IV. la) 



A(r,t) 



27r X - 1 r ^ 

* dt' 



iV 



i(cqt— qr) 



^||dVe qj .[J P ](r',t')e^' + -') 



+ e qiJ e* ( - c " t+q - r) / ^ / dVe^ ■ [J P ](r', t'y^' 



then we can re-write Equations (jllLlll Hll.91 and IIH.7a)l as follows: 



v(t) = / d 3 r 



p T (r)0(r,t) - -J T (r,t)-A(r,t) 

c 



It follows immediately from these definitions of <p(r,t) and A(r,t) that 

V 2 ^(r,t) = -4vrp P (r,t) 
V-A(r,t) = 0. 



(IV. lb) 



(IV.2) 



(IV.3a) 
(IV.3b) 
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The latter equation is a consequence of our use of transverse photon modes (e qj - • q = = 

To proceed, we replace the sum by the integral ^ f d 3 q. Moreover the completeness 
of the set (e q ,i, e qj2 , q) allows us to write 

2 

J P = q (q • J P ) + ^ Kd ' J p) 

from which we can obtain 

q(q • J P ) 



+ < '<-■"-<*■" / df / dV ( J P (v',t') - q(q- J P( r, »*0) ) ^(oBf-q-rO 



5^ ejj (eqj ■ J p) = % ■ Jp) = Jp - q (q • Jp) = Jp - 

J=l 3=1 

Then we can rewrite Equation ()IV.lb|) as 
A(r,t) = JL y £i [ e ^') jf" ^ df I dV (jp(r',i') - q(q- y /;t/)) ) e ^(-^'+ q -') 

(IV.4) 

y- / j 

By carrying out the differentiations, one can show that 

J" (*™ - ) — 

(IV.5a) 

and 

But projectile charge conservation implies that 

= J dVJ P (r',t) • V r /e iq '( r - r ') = -i J rf'Vq- J P (r , ,t)e iq - (r - r ' ) . 
Thus Equation (|IV.5bJ) becomes 

|v«r, .) = JL / A / d V ^-y^W >, (,V.5c) 
which we can combine with Equation (|IV.5a|) to get 

V»A(r,«) - - i|v«,,«) = - JL / dV W «) / ») 
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= -«V / rfV (27r) 3 5(r-rOJp(r',t) 

V 2 A(r, t ) - p -gLlJ _ -_v v( r, t ) = --J P (r, *) (IV.6) 

Equations (|IV.3al IIV.3bl and IIV.6|) show that A(r, t) and <^(r, t) defined in Equations 
(jlV.lal IIV.lb|) are the classical vector and scalar potentials A(r, t) and <p(r,t) associ- 
ated (in Coulomb gauge) with projectile charge and current densities (pp(r, t), Jp(r, t)) . 
Morevover, Equation (jIV.2)) shows that the effective potential v(t) is the classical interac- 
tion of ((p(r, t), A(r, t)) with the target charge and current densities. 

Thus we have proven that an RCE calculation in which the classical electromagnetic field 
of the projectile is used to induce transitions between target states gives the same target 
transition amplitudes, to all orders of perturbation theory, as would a calculation in which 
the interaction between projectile and target is mediated by a quantized electromagnetic 
field. 



V. THE LORENTZ GAUGE 



The two previous sections used the Coulomb gauge. In this section, we will describe the 
modifications needed if the Lorentz gauge is used. 

1. In the Lorentz gauge, the full Hamiltonian is 

H = h + J d?r(^ (p P (r, t) + p T (r)) y?(r) - i ( J P (r, t) + J T (r) ) • A(r) ) , (V.l) 
which lacks the density-density interaction present in Equation (jIII.5aj) . 

2. In the Coulomb gauge, the vector potential A has only two components, and they 
are transverse. In the Lorentz gauge, A has these two tranverse components, and 
also a longitudinal component (along q). Moreover, the scalar potential ip is also 
quantized. Thus, in addition to the two transverse photons of the Coulomb gauge, 
we have a longitudinal photon and a scalar photon. The commutation relations of 
the longitudinal photon creation and annihilation operators are the same as for the 
transverse photons (Equations (jIII.3a|) and (|III.3b|) ). However, the treatment of the 
scalar photon is more complicated. A consistent formalism for quantizing the scalar 
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field was developed by S.N. Gupta 



and K. BleulerP]. For our purposes, the only 
manifestation of the extra complications of scalar field quantization is the presence of 
the two minus signs in Equation (V.2) below. 



When these changes are accounted for, Equation (jIII.9|) is replaced by 

q 



(V.2) 



j—l J J tp J 



e^ r 3 P (v',t' p )e^ c < + ^ 



d 3 re^ • [J T ] MI/ (r)e^-^ +t »- r ) J J d 3 r% d • J P (r', t' p )e^ c <-^ 



- c 2 / d 3 r[p T }^(r)e 



r°° tit' r 

rf 3 r[p T ]^(r)e i ^ c ^ +q - r ) J J d 3 r'p P {r', Qe^'p'^ 

Following the procedure of Section 4, we cast this expression in the form (jlV.2|) . which leads 
to 

3 



<p(r,t) 



2C7C 



y-y-Ir^-q-r) f°° dt , f rf 3 r / pp / r / )t / )e i(- Cg t'+q-r') 

t 



A(r,t) = 



+ e i( - c<? * +qr) / dt' I rf 3 rXr / ,t , )e i(c9 *'" q "' ) 
2tt 



(V.3a) 



q j = l q Jt J 

+ e qii e i( - C9t+q - r) / dt' I d 3 r't^ ■ [J P ](r', *yM'-<i- r ') 
The polarization sum simplifies to 

3 



{-cqt'+q-r') 



(V.3b) 



and the replacement of ^ by the integral ^ J d 3 q still applies. When these relations are 
used, it is staightforward to calculate that <p(v,t) and A(r, t) defined by Equations (jV.3al 
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and IV. 3b]) satisfy 



V 2 



c 2 dt 2 



^ 1 d 2 
V 

c 2 dt 2 



<f( r ,t) = -47rp P (r,t) 
A(r,t) = J P (r,t) 



V- A(r,t) + 



l d<p(r,t) 
c dt 







(V.4a) 
(V.4b) 
(V.4c) 



Equations ()V.4al IVT4p1 and lV.4c|) show that <p(r, t) and A(r, t) defined in Equations ()V.3al 
IV.3b|) are the classical vector and scalar potentials A(r, t) and (p(r, t) associated (in Lorentz 
gauge) with projectile charge and current densities (pp(r, t), Jp(r, t)) . Thus, whether one 
chooses to work in the Coulomb or Lorentz gauge, in an RCE calculation of target transition 
amplitudes the effect of the quantized electromagnetic field can be replaced by the classical 
electromagnetic field of the projectile used in the classical expression IIV. 2l for the interaction. 

The proof of the equivalence of the quantized field and classical field treatments of 



Coulomb excitation given in ABHMW |4| applies only to the on-shell (u = oo ja ) Fourier 
component of v ia (t). The reason for this apparent restriction is that ABHMW used the 
Coulomb gauge for the calculation of the quantized- field version of v 10l {t), but they used 
;he Lorentz gauge for the calculation of the classical field version. It is shown in Reference 
1^ 1 that the versions of v ia (t) calculated in the two gauges agree in their on-shell Fourier 
components, but are generally different when uj 7^ u^ Q . In Sections III, IV, and V, we have 
shown that the quantized field and classical field treatments agree for all t, and so for all uj, 
if both are calculated in the same gauge. 



VI. DISCUSSION 



It has been emphasized that an important ingredient in our derivation is the assumption 
that the charge and current densities associated with the projectile are specified functions of 
position and time, which is equivalent to the assumption that the projectile does not change 
its internal state during the collision. This assumption is also made in RCE calculations 
involving classical electromagnetic fields. Even if both the target and projectile are allowed 
to undergo internal transitions, the fields that induce those transitions are always assumed to 
be generated by moving static spherically symmetric charge distributions, with no internal 
structure. To remove this assumption, we would have to calculate, at each instant, the 
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classical retarded electromagnetic fields that have been generated by the projectile and 
target since the beginning of the collision, while they have moved relative to each other and 
undergone changes of their internal states. These would be the electromagnetic fields that, 
in turn, would be used to induce further internal transitions. It is clear that this would be 
a very difficult calculation, even though it involved only classical fields. It would be the 
classical analogue of the quantized field calculation we have described in Sections II and III, 
if we had allowed the projectile state to change as well as the target state. 

Although it would be difficult to perform the complete calculation just described, in 
which the electromagnetic fields are generated by the actual dynamic transition charge and 
current densities, it may be required for a full understanding of ultra-high-energy RCE. 
Baltz, Rhoades-Brown and Weneser j^] have estimated excitation probabilites for collisions 
between oppositely directed 100 GeV Au nuclei. They find excitation probabilities that are 
greater than 0.5 for grazing collisons, and greater than 0.1 out to impact parameters of about 
50 fm. Thus it is important to investigate, with classical and/or quantized electromagntic 
fields, the effect on excitation probabilities of internal transtions occurring within both the 
target and projectile. 
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Figure Captions 

Figure 1. The ^0(^2)^0(^1) term of Equation (|III.12|) . Time increases as we move upwards 
in the diagram. The vertical line on the left corresponds to the projectile, and the lines 
labelled a, j3, 7 represent states of the target. 

Figure 2. The w(ti)vo(t2) and Va(ti)w(t2) terms in Equation (|III.12|) . The curly lines 
represent exchanged photons. The t[ and t' 2 variables must be integrated from —00 to 00. 

Figure 3. The wftijwfa) term of Equation (jIII.12|) . The t[ and t' 2 variables must be 
integrated, independently, from —00 to 00. 



17 




Fig .1 




Fig .3 



